Entanglement purification protocols for all graph states 
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We present multiparty entanglement purification protocols that are capable of purifying arbitrary 
graph states directly. We develop recurrence and breeding protocols and compare our methods 
with strategies based on bipartite entanglement purification in static and communication scenarios. 
We find that direct multiparty purification is of advantage with respect to achievable yields and 
minimal required fidelity in static scenarios, and with respect to obtainable fidelity in the case of 
noisy operations in both scenarios. 
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I. INTRODUCTION 

Entanglement purification (distillation) is an impor- 
tant primitive in quantum information processing. It al- 
lows one to overcome the influence of noise in quantum 
communication P, 0, 0] , and enables one to obtain prov- 
able secure quantum key distribution even in the con- 
text of noisy channels and over arbitrary distances 
0. The applications of entanglement purification are 
however not limited to bipartite communication scenar- 
ios. First, it has been extended to certain multiparty 
scenarios [1, 0, H, S leading e.g. to novel quantum 
primitives such as multiparty secure state distribution 

Second, in the context of fault-tolerant quantum 
computation, entanglement purification is a key ingre- 
dient to obtain improved error thresholds [H, [ll. Ap- 
plications in quantum error correction (T^ and quantum 
simulation fl5l| have also been discussed. 

The basic idea of entanglement purification is to use 
several copies of a noisy entangled state to generate, by 
means of local operations and classical communication 
(LOCC), a few copies with improved fidelity. So far, 
entanglement purification protocols have been developed 
that are capable of purifying Bell states [l|, 0, 0] in the 
bipartite case, and all two-colorable graph states d, Q 
as well as W states 16] in the multiparty case. They 
also have been generalized to higher dimensional systems 

Graph states [13, are a family of multiparty 
entangled states with interesting entanglement proper- 
ties. Two-colorable graph states are sub-family associ- 
ated with two-colorable graphs and include a number of 
interesting states such as GHZ states, cluster states and 
codewords of CSS error correction codes. Graph states 
appear for instance in the context of measurement-based 
quantum computation, where a given graph state rep- 
resents an algorithmic-specific resource that allows one 
to realize a specific unitary operation on several qubits 
by local measurements. Typically, these graphs are not 
two-colorable, but one may wish to purify this resource, 
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e.g. to realize one-way quantum computation in a fault- 
tolerant manner. 

In this paper, we present entanglement purification 
protocols (EPP) that are capable of purifying arbitrary 
graph states. To be precise, we develop for each graph 
state a direct multiparty recurrence protocol and a mul- 
tiparty breeding protocol. Our key ideas are twofold. 
First, we show that if auxiliary (even noisy) two- colorable 
graph states are available, the fidelity of the noisy graph 
states can be improved. Second, we show how to get 
a single copy of the auxiliary two-colorable graph state 
from two identical copies of the fc-colorable graph state 
by LOCC. Thus, unlike the known two-colorable graph 
state entanglement purification protocol [1, toj , our pro- 
tocol utilizes different shapes of graphs. 

The protocols are applicable both in (i) a static 
(LOCC) scenario, where the parties attempt to purify by 
LOCC several copies of given noisy multiparty entangled 
states, as well as in (ii) a communication scenario where 
the parties are allowed to generate (arbitrary) multiparty 
states locally, distribute them through noisy quantum 
channels and attempt to end up with high-fidelity tar- 
get graph states. Our first idea is commonly utilized to 
both scenarios. The second idea is crucial in the static 
(LOCC) scenario, while in the communication scenario 
we may prepare separately auxiliary two-colorable graph 
states through noisy channels in an effective way. 

In both cases, we show that the new direct multiparty 
EPP are superior to alternative approaches based on bi- 
partite entanglement purification. In particular, we find 
that in the static scenario (i) the yield of direct multi- 
party breeding protocols is higher than for any strategy 
based on bipartite purification, and also the purification 
regime is larger. When also considering noisy local con- 
trol operations, we show that for both scenarios (i) and 
(ii), the new multiparty recurrence protocol allows one 
to reach higher fidelities. We first review the concept of 
graph states in Sec. Ill Al and then present new recur- 
rence and breeding protocols in Sec. IIV Al and [TVBl We 
present an alternative strategy based on purification of 
two-colorable sub graph-states in Sec. |Vl A comparison 
with bipartite strategies for noiseless in Sec. IVI Al and 
noisy local control operations in Sec. lVIBl finallv demon- 
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strates the advantage of these new protocols. 

II. GRAPH STATES AND MANIPULATION 

In this Section, we summarize the basics concerning 
graph states and their manipulations. 

A. Definition and notation for graph states 

Graph states are a family of multiparty entangled 
states associated with mathematical graphs fli*!. A 
graph G = {V, E) is given by a set ^ = {1, 2, ... , N} 
of N vertices connected in a specific way by edges E. 
To every such graph there corresponds a basis of N- 
qubit states {\h)g}, where each of the basis states 
{ix = /LtiyU2 . . . /kat) is the common eigenstate of N com- 
muting correlation operators with eigenvalues (—1)''" 
such that = 0, 1. That is, they fulfill the set of eigen- 
value equations K'^\^)g = (— 1)'^° |/x)g, a — 1,...,N. 
The correlation operators are uniquely determined by the 
graph G and are given by 

{a,b}GE 

where cri"'' denotes the application of the corresponding 
Pauli operator (a — x, y, z) by the party a. Equivalently, 
any graph state \ij.)g can be written in the following man- 
ner: 

iM)G=n(^i'^'r" f n Az'"^)] (2) 

a=l \{b,c}eB / 

where KZ — diag(l, 1,1,-1) in the computational basis 
is the controUed-phase gate, and |±) — -^(|0) ± |1)). 

We will use the concept of fc-coloration in the follow- 
ing. A graph is called fc-colorable if there exist k sets of 
vertices Ai , A2 , . . . , e V such that there are no edges 
within each of the groups Aj for all j, i.e. for all a, 5 G Aj , 
and for all j, we have {a, 6} ^ E. Two-colorable graphs 
are a special instance with k — 2, where multiparty en- 
tanglement purification protocols are known @,l9| . How- 
ever, only a subset of graphs is two-colorable and, in 
principle, a graph may be A^-colorable. We remark that 
local-unitary equivalent graph states [ll| may correspond 
to graphs with different coloring, and the minimum k 
within the local equivalence class is not known. Under 
local Clifford operation, we have however that generally 
fc > 2, i.e. not all graph states are locally equivalent to 
two-colorable graphs. 

Associated with a given fc-colorable graph G with col- 
oring {Ai,A2,...,Ak}, we define two-colorable graphs 
{51, 52, ■ ■ • , (see Fig H] for an illustration), where gj 
contains only the edges between the set Aj and the re- 
maining sets {Ai,i j}, but where edges between the 
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FIG. 1: (Color online). A 3-colorable graph G and the 3 
corresponding two-colorable sub-graphs gi, (72, and 513. gi 
corresponds to the red color (vertices 2,4 and 6), g2 to the 
green color (vertex 7) and gs to the blue (vertices 1,3 and 5). 

remaining sets are erased. That is, the sets {Ai,i ^ j} 
form a new set A^ = V\Aj. The set of indices fj,a corre- 
sponding to the set Aj will be denoted by fij. Note that 

U,ti<?, = G. 

In this paper, we consider mixed states diagonal in the 
graph-state basis, 

PG^Yl ViM2...Mfc IM1M2 ■ • ■ Mfe)G(MlM2 ■■■^J■k\, (3) 

where we have grouped the multi-index /x into k multi- 
indices fij — /ijj^ . . . fj,j^ (m = \Aj\) corresponding to the 
sets Aj defined by a chosen fc-coloration of the graph G, 
note that mixed states resulting from any noise models 
can be brought to this form by means of local depolar- 
ization, i.e. by applying randomly the local operations 
corresponding to the correlation operators {K^}. Diago- 
nal elements in the graph-state basis, i.e. the coefficients 
A^i/X2 .-Aifc remain unchanged by this procedure, and any 
mixed state can be assumed to be of the form ^ without 
loss of generality. Hence we can interpret the mixed state 
as an ensemble of graph states |/.t)G, where \iji)g appears 
with probability A^. We can therefore restrict our at- 
tention to the transfer of these indices between unknown 
pure states as presented in the next subsection. 



B. Operations on graph states 

We briefly mention two operations which play key roles 
in entanglement purification protocols. One is the multi- 
lateral CNOT which enables one to transfer the stabilizer 
eigenvalues fi between two states. The other are Pauli 
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FIG. 2: (Color online). Effect of a CNOT gate on a product 
of two graph-states. The final state is associated to a graph 
composed of the two initial graphs with additional edges be- 
tween the control qubit 61 and all neighbors of the target 
qubit ti2. 



measurements which allow one to evaluate /x with the 
help of classical communication. 

Let us first describe the action of a CNOT gate on a 
product of two (possibly different) A^-qubit graph-states 
\h)gi and \v)g2- We apply the CNOT gate from ai (a*'' 
qubit of the first state) to 02 (a*'' qubit of the second 
state). We deduce the expression for the resulting state 
from ([2]). As |-|-) (gi |-|-) is an eigenvector of CNOT asso- 
ciated to the eigenvalue 1, the difference between the ini- 
tial and the final state is due to the commutation relation 
between the CNOT gate, the cr^ operators and the con- 
trolled phase gates. A straightforward calculation yields 



CNOT^'^i^'^^) |/x)gJiv)g^ = 



(4) 



where Na^ are the neighbors of vertex 02. The final state 
is related to a graph of 2 A'^ vertices which is composed of 
the two initial graphs with addition of all edges between 
vertex ai and the neighbors of vertex 02. See Fig[2]for an 
illustration. In addition, bit is fiipped if bit = li 
which corresponds to a transfer of information from the 
second state to the first one. Note however that this par- 
ity information cannot be evaluated without disturbing 
the first state. 

To read out the information transfered between two 
states, one uses Pauli measurements in directions that 
correspond to the correlation operators {K^}- Suppose 
we want to determine the bit string fjbj corresponding to 
a color Aj for a given graph state \fj,j,^j)Q. Then, the 
parties belonging to color Aj measure their qubit a in 
the eigenbasis {|±) = :^(|0) ± |1))} of a^, obtaining re- 
sults G {0, 1}, while all the other parties make their 
measurement in the eigenbasis {|0), |1)} of Cz, obtaining 
results (b G {0, 1}. According to Eq. dT]), each stabilizer 
eigenvalue fia of {a G Aj) can be evaluated via clas- 
sical communication of these measurement results as. 



a — £.a Cb 



(5) 



Note that all indices corresponding to color Aj are mea- 
sured simultaneously. 



III. STATIC SCENARIO AND 
COMMUNICATION SCENARIO 

In this article, we are interested in the generation of 
high-fidelity graph states shared among N spatially sep- 
arated parties. We consider two cases (i) a static LOCC 
scenario, and (ii) a communication scenario, since both 
situations are of practical relevance in quantum informa- 
tion processing. 

In (i), the static LOCC scenario, the parties share M 
copies of a mixed state p. They can manipulate the states 
by means of LOCC in order to generate a high fidelity ap- 
proximation of the target state \0)g- No quantum com- 
munication between the parties is allowed. The efficiency 
of this procedure is measured by the yield, i.e. the ra- 
tio of high fidelity output state per input state. In (ii), 
the communication scenario, the parties are connected 
by noisy quantum channels and are allowed to distribute 
any locally generated entangled state through the noisy 
quantum channels, in addition to sequences of LOCC. 
The efficiency of this procedure is measured by the (in- 
verse of the) quantum communication cost |19|, i.e. the 
number of channel usages required to generate a desired 
multiparty entangled state with sufficiently high fidelity. 

The scenario (i) deals with entanglement properties of 
a given mixed state p, i.e. whether high fidelity multi- 
party entangled states can be distilled from several copies 
of p by LOCC. The scenario (ii), on the other hand, has 
a certain physical set-up (a communication set-up) in 
mind, and hence deals with the question whether in such 
a context high-fidelity entangled states can be generated. 
For our purposes, the most relevant difference between 
the two scenarios is that in case of (ii) noisy entangled 
states of any kind can be generated by distributing them 
through noisy quantum channels, which can be used to 
generate the desired target state. In (i), only copies of 
p are given which should be manipulated by LOCC. If 
different states are needed, this has to be prepared by 
manipulating several copies of p by LOCC. 



IV. DIRECT MULTIPARTITE 
ENTANGLEMENT PURIFICATION PROTOCOL 

A. Recurrence protocol 

We now present a recurrence protocol that allows one 
to purify directly any A;-colorable graph-state, provided 
auxiliary two-colorable graph-states corresponding to the 
k different colors are available. The sub-protocol is ap- 
plicable in both scenarios, and we describe later on how 
to construct the auxiliary two-colorable graph-states in 
each scenario. Note that if the target graph state is two- 
colorable, auxiliary graphs 171 and 32 are nothing but 
G. Thus our protocol covers known protocols for two- 
colorable graph states. The total protocol consists of k 
sub- protocols {Vj} {j = 1, . . . , fc), each of which serves 
to purify a state partially, i.e. with respect to the index 
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vector fij corresponding to a group Aj (i.e. a certain 
color). We outline the protocol in the following. 

1. Sub-protocol Vj 

Assume that auxiliary (possibly noisy) two-colorable 
graph-states pg- are available in addition to an ensemble 
of fc-colorable graph-states pc- We purify the fc-colorable 
graph-states pc with respect to bit string jij: 

(1-1) The parties take two noisy states in such a way that 
the first state is po and the second one is pg- . The 
parties in the group Aj apply CNOT gates from 
the second (control) state pg. to the first (target) 
state PG, while the parties in the group A^ apply 
the CNOT gate in the opposite direction. 

(1-2) The parties measure the second state Pg. locally in 
the eigenbasis of Ux for the group Aj and in the 
eigenbasis of (j^ for the group Aj. Using classi- 
cal communication, they decide to keep the first 
state Pg if all eigenvalues corresponding to the cor- 
relation operators {K^,a e Aj} of the group Aj 
(determined by Eq. ([5])) are 0, or to discard it oth- 
erwise. 

Let us take a close look at the sub-protocol Vj. In (1-1), 
the multilateral CNOT operation, 

]J CNOT^'^^^'^i) J]^ CNOT^''!^''^), (6) 

{ai,a2}eAj {bi,&2}eAj 

is applied between a first state, corresponding to the 
graph G, and a second state corresponding to the graph 
gj (from now on, we use ai for the vertex of the l-ih 
state in the group Aj and hi for the one in the group 
Aj). A straightforward calculation using Eq. (jH) gives 
the following map 

\fJ.j,flj)G\Vj,Vj)g, ^ \flj,fJ,j ®Uj)GWj ® fJ-j,V-j)g,, (7) 

where © means bitwise addition modulo 2, which shows 
the transfer of information about the stabilizer eigenval- 
ues between the two states. Note that the final states 
correspond to the same graphs as the input states. After 
that, one measures the second state in order to determine 
the indices i/j © fij corresponding to color Aj , using the 
procedure described in (1-2), where each bit i^j (Bpj is de- 
termined using Eq. ((S]). If all the parities Uj^^j are 0, it 
is expected that pa — ^ and t'a = for a G Aj are prob- 
able since |0)g and |0)g^. have been assumed to be the 
majority in their ensembles. That is why then the first 
state is kept and otherwise discarded. As consequence, in 
the expansion ([3]) of the density matrix, elements of the 
form Ao,-y~. are increased. One finds that the new matrix 
elements of p'q are given by 



where k is a normalization constant guaranteeing that 
tr(p') — 1, and A are the coefhcients for the two-colorable 
graph-state Pg^ written in the form ([3]) . We remark that 
here we do not address the question of the unfavorable 
scaling behavior of the efficiency of the proposed protocol 
with the total number of particles N, for recent devel- 
opments on this subject in the context of two-colorable 
graph state purification see [20| . 

We illustrate how the protocol works by looking at the 
simple toy-case where noise acts only on one color. Let 
us consider G as the 5-qubit ring and gi as the 5-qubit 
cluster, and study the effect of Vi on mixed states of the 
form 

PG = J2 Vi AO lAti, 0, 0)G(Ati, 0, 0|, (9) 

Ml 

Pgi =^VAo|i^i,O,0)g,(i/i,O,0|. (10) 

Note that even though gi is two-colorable, we group the 
vertices in three distinct sets corresponding to the col- 
ors of the ring. As only the qubits in Ai are noisy, 
a sequence of applications of sub-protocol Pi is suffi- 
cient to purify the state. A straightforward calcula- 
tion gives the new coefficient of the purified ring-state: 

•^71,0,0 ^ A-y-^^o.o A-y.j^,o,o/ X^i/i (^j'l.o.o Ai.,-^,o,o)) meaning 
that the dominant coefficients are increased. 

The whole purification protocol consists of a sequence 
of applications of the sub-protocols Vj corresponding to 
all colors j — 1, . . . ,k. Even though there is a back-action 
of noise for the colors which are not purified for the step 
j, one obtains an overall increase of the fidelity Ao if the 
fidelity of the initial state is sufficiently high. In fact, 
Ao = 1 is an attractive fix point of the protocol under 
the ideal local operations. 



2. Preparation of the auxiliary state pg. in static scenario 

Next we describe how to obtain the auxiliary two- 
colorable graph-states pg. required for sub-protocol Vj. 
In the communication scenario, these states can be gener- 
ated directly by distributing them through (noisy) quan- 
tum channels. In the static scenario, the situation is 
slightly more complicated and a pre-processing to the 
purification sub-protocol Vj is needed which we will de- 
scribe in the following. 

Assume that an ensemble of noisy fc-colorable graph 
states PG is available. We supply an ensemble of the 
auxiliary two-colorable graph state pg. by LOCC: 

(0-1) The parties take two identical noisy copies of pc- 
The parties in the group Aj apply a CNOT gate 
from the second (control) copy to the first (target) 
copy, the parties in the group Aj apply the CNOT 
in the opposite direction. 
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(0-2) The parties measure the second state in the eigen- 
basis of CT2. By this, they erase all the edges be- 
tween the two states and they are left with a state 
corresponding to the two-colorable graph gj after 
a suitable change of local basis depending on the 
measurement outcome. 

A straightforward calculation using Eq. ^ shows that 
the multilateral CNOT applied in step (0-1) results in 
the following map when applied to a product of two k- 
colorable graph-states 

^ II AZ(''''''-^\fij,fij(Buj)g^\uj®fij,u-)G, (11) 

which is a product of a graph-state associated to gj (state 
1) and a graph-state associated to G (state 2), with addi- 
tional edges between the two graphs. The local cr^ mea- 
surements on the second state (step (0-2)) erase the cor- 
responding vertices and all edges associated to them in 
the graph. This ensures that after the measurement the 
first state is a two-colorable graph state corresponding to 
the graph gj as desired. 

Note that the multilateral CNOT does not only create 
the desired two-colorable graph-state, it also results in a 
transfer of information between the two states. In par- 
ticular, the part of the index bit of state 2 corresponding 
to color Aj is given by Vj © fXj, where and i' corre- 
spond to the states before the CNOT operation. We are 
therefore not only able to create the two-colorable graph- 
state but also to perform at the same time a first step of 
purification. To this aim, we replace the measurement in 
the eigenbasis of ct^ by a measurement of the correlation 
operators with a & Aj, where we keep the state only 
if the expectation values of all these correlation operators 
are zero. That is, step (0-2) can be replaced by 

(0-2)' The parties measure the second state locally in the 
eigenbasis of ax for the group Aj and in the eigen- 
basis of CTz for the group Aj . The first state is kept 
if all eigenvalues corresponding to correlation oper- 
ators {K^,a G Aj} of the group Aj (determined 
by Eq. ©) are 0, and discarded otherwise. 

Note that depending on the measurement outcomes a lo- 
cal unitary operation on the remaining copy is required 
to ensure that \0)g- is the dominant component of the 
resulting state Pg^ . To be precise, the parties should ap- 
ply a local unitary operation ObieAj, b'^eN^^nA-i^^^'^)''''''' 
depending on their measurement pattern (^1,1 after either 
(0-2) or (0-2)'. 

B. Breeding and hashing protocols 

Hashing and breeding protocols were introduced for 
the bipartite Bell state in Ref. [l|, 0| , for the GHZ state in 



Ref. [3, and generalized to all two-colorable graph states 
in Ref. E i. E HII • We will now show that the muhilat- 
eral CNOT operation together with the use of different 
states allow one to construct hashing and breeding pro- 
tocols for any graph state. In both cases, we are given an 
ensemble of M imperfect fc-colorable graph-states, with 
M — » 00. One then transfers information from randomly 
chosen subsets of rn states to perfect states (breeding) or 
imperfect states (hashing), which are then measured. At 
each round our knowledge about the remaining states is 
increased. 

We propose here a generalization of the breeding pro- 
tocol to all graph-states. A hashing protocol can be ob- 
tained in a similar manner with the additional require- 
ment of taking into account the back-action due to the 
imperfection of the states used to read out the informa- 
tion. One is given M copies of an A^-qubit graph-state 
corresponding to graph G, where G is fc-colorable. We 
consider M — > 00. In addition, one is also given k en- 
sembles Ej, j — 1, . . . , fc, which one needs to give back 
at the end, where ensemble Ej contains perfect copies 
of the two-colorable graph-state corresponding to graph 
gj . Let us call Bi, i £ 1 , . . . , the binary vector which 
contains the value of bit i for a sequence of m states. 
It is possible to determine the parity of all Bi belonging 
to a given colour Aj simultaneously by measuring only 
one state. This is done by performing CNOT gates from 
one state of Ej to the m fc-colorable graph states for all 
qubits in Aj and in the opposite direction for the other 
qubits. To recover the parity of B.^, one measures the 
two-colorable graph-state using the method described in 
Seq. IIIBI By this, one determines the eigenvalues of all 
correlation operators corresponding to set Aj, which are 
given by lU . One needs to repeat the procedure at most 
M S{a[^\a[^'') times to obtain perfect knowledge of Bi, 
where S{af'\a[^^) = — a-*^^ logj a-°^ — logj a-^^ is the 
entropy and a[^''' = E^^^^^. A^,^,...^,...^„ . 

For each of the fc colors, one performs the sequence 
of operations described above ma,Xi^Aj{M S{a[^\a^^^)) 
times in order to obtain the parity of all Bi belonging to 
this color. Given this information, one ends up with a 
pure state corresponding (up to local unitary operations) 
to M' copies of |0)g. The last step consists in reconstitut- 
ing the pool of perfect two-colorable graph-states we were 
given at the beginning. This is done by performing the 
multilateral CNOT given by ^ on two copies, using by 

this 2 maxig^^. (Af S{af\a[^^)) fc-colorable graph-states 
for each color. A lower bound Y of the yield is therefore 
given by 

fc 

y = 1 - 2 ^ max.eA, S{a,). (12) 
i=i 

We remark that this construction might not be opti- 
mal, and a significant improvement of the yield could 
be achieved by an optimal procedure that generates M' 
copies I'^gj) from M copies of |^'g)- To illustrate our 
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FIG. 3: Yield as function of fidelity for the breeding protocol 
applied to the 5-qubit ring-state. 



protocol, we calculated the yield for a 5-qubit ring state 
of the form 

p = f \0)g{0\ + (1 - /)/(25 - 1) (1 - |0)g(0|) . (13) 

Note that this state is a 3-colorable graph-state. The 
entropy S{af'\af'^) is identical for all bits. It is given by 

-(/+(2^-l)^)log.(/ + (2^-l)^) 

The yield, given by (1 — 6 S{a^^^ , aj^^ )) is plotted as func- 
tion of / in Fig.[3l One sees that the yield is approaching 
1 for a state of fidelity close to 1. 

V. ALTERNATIVE PURIFICATION 
PROTOCOLS VIA CUTTING INTO AND 
RECONNECTING 2-COLORABLE GRAPH 
STATES 

Until now we have considered strategies where the k- 
colorable graph-state was directly purified. Another pos- 
sibility is to purify smaller parts of the graph-state corre- 
sponding to two-colorable sub-graphs and connect them 
at the end. The advantage of this strategy is to make pos- 
sible the usage of known two-colorable graph state purifi- 
cation protocols. Different intermediate strategies can be 
designed depending on the chosen sub-graphs. We use an 
iterative method to construct a set of sub-graphs with no 
overlapping edges, such that their union is the graph G. 
Step j of the procedure consists of choosing color Aj out 
of the k ~ j remaining colors and constructing sub-graph 
gj, such that gj is a further sub-graph gj of Sec. Ill Al 
where all vertices a G Ai,i = 1, . . . , j — 1 have been 
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FIG. 4: (Color online). A 3-colorable graph G with colors 
Al = {7}, A2 = {2,4,6} and A3 = {1,3,5}. 51 and 32 
are two-colorable sub-graph which can be obtained from G 
and which give G when merging them together. We use the 
method described in Sec.|V]to construct the sub graphs. For 
example, to obtain gi, we choose Ai as first color and erase 
all edges between A2 and A3 using the procedure described 
in Sec. IIV A 21 The sub-graph (72 is obtained by erasing the 
vertex Ai via a az measurement, and no further processing is 
required since this graph is already 2-colorable. 

erased. At each step, one checks if the graph G\{Ul^j^gi} 
is two-colorable. If this is the case, one stops the proce- 
dure. The states corresponding to the sub-graphs are 
distributed from the beginning in a communication sce- 
nario, while they are generated from two copies using the 
multilateral CNOT operation described in Sec. IIV Al in a 
static scenario, followed by a measurement in the az basis 
of the vertices one wants to erase. Once the states are dis- 
tributed, the multipartite purification protocol for two- 
colorable graph-states [1, 0| is applied. To conclude the 
protocol, one merges the different graph states together 
in order to create the final fc-colorable graph-state. Given 
two vertices ai and 02 belonging to different graphs, the 
corresponding party merges them by applying a projec- 
tive measurement given by Pq = |0)(00| -I- |1)(11| and 
Pi = |0)(01| + |1)(10| (with outcomes and 1). A local 
correction given by Jlbe-ZV '^^^'^ applied to the state 
when the measurement result is 1. 



VI. COMPARISON WITH BIPARTITE 
STRATEGIES 

A. Performance under ideal local operations in 
static scenario 

1. Yield 

We show in this section that in a static scenario, di- 
rect multipartite purification is more efficient than strate- 
gies based on bipartite purification. This is due to the 
fact that any strategy based on bipartite entanglement 
purification requires that at a certain stage Bell-pairs 
shared among pairs of parties are generated. After the 
entanglement purification protocol, another sequence of 
LOCC must be applied to recover the multiparty entan- 
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gled states. We show here that these two sequences of 
operations necessarily generate losses even when applied 
to pure states. We illustrate this by considering as exam- 
ple the 5-qubit ring state to which we apply the method 
introduced in Rcf. 9] to quantify the loss. Wc start with 
an ensemble of M perfect ring states |0)g(0|, which are 
then transformed to an ensemble of Bell pairs |0)''g^\o|, 
shared between the different parties, by means of LOCC. 
Another sequence of LOCC brings the pairs to an en- 
semble of M ring states. The total procedure can be 
summarized as follows: 

|0)g(0|»*'^ _ (g) |o)(g(0|»'"- ^ |0)g(0|«^^. (15) 

(kd)-M<l 

We now calculate a bound on the yield M jM for this 
procedure. To do it, we apply the following inequalities 
which were used in Ref. [23] to show the irreversibility of 
entanglement transformation between singlets and GHZ 
states: (a) The entropy can only decrease on average 
under LOCC operations; (b) The average increase of rel- 
ative entropy of system B = V \A\s smaller than the 
average decrease of entropy in system A. If we consider 
a density operator pab describing a pure state which is 
transformed to an ensemble {pfc,^^^} by LOCC opera- 
tions we have that (a) 

5(/)>^p.5(p^), (16) 
fc 

where S{p"^) is the von Neumann entropy of the reduced 
density operator for system A, p is the initial state and p 
the final state. We consider the bipartition of the system 
of qubits into system A and system B — V \ A. Then, 
we have (b) 

Y^pMpf) ~ E^ip^) < Sip^) - 5]p.5(p^), (17) 

k k 

where Er{p^) is the relative entropy of entanglement for 
p^ 

Erip^)^ min 5(/||a^), (18) 

fT^scparablc 

with 

S{p^\\a^) := Tr(/log/) -Tr(/loga^) (19) 

being the relative entropy. Let us now apply this in- 
equalities to our example. We first calculate (fT6|) for 
the second part of process with A = {01,02} and 
B = V \ A. We have for the state describing the ensem- 
ble of pairs Sip"^) = J2b:beB ("^aif* +'ma2b)- This comes 
from the fact that the entropy of a pure state is zero and 
that Ttb p is not a pure state if and only if a single of the 
2 qubits belonging to set A is part of an entangled pair. 
In this case, the entropy is equal to unity for a single 
state. As there are mab states containing a fully entan- 
gled pair between qubits a and b, the contribution of this 



states to the total entropy is niab- For the M ring states 
we find S{pa) = 2M, as TrB(|0)G(Ol) can be written as 
a sum of 4 projectors with equal weights (see below), and 
hence for a single copy of the ring state, the entropy of 
entanglement with respect to the bipartition in question 
is two. Summing up the contributions of all bipartitions 
of 2 and 3 qubits, we obtain 

6 ^ mab > 20 M. (20) 

a,6;a<fc 

We now apply inequality (fT7|) to the first part of pro- 
cess ([T5t . with A = {01,02} and B = V \ A, to 
have a bound on M. We distinguish two kinds of 
bipartitions: (i) qubits oi and 02 are neighbors (ii) 
they are not neighbors. Let us study the entangle- 
ment properties of the state |0)g(0| obtained by tracing 
out A. In case (i), Tr^(|0)G(0|) = 1/4 (|000}g'(000| + 
|001)g'(001| + |100)g'(100| + |101)g'(101|), where G" 
stands for the graph corresponding to the 3-qubit GHZ 
state. A straightforward calculation shows that this 
state is separable as it can be written as Tr^(|0)G(0|) = 
(Had(g)l ® Had)''' p' (Had®l ® Had), with p' = 1/4 (| + 
-l-f)(+ + +1 + I - -+)(- - +1 + I - +-)(- + 

— I + I + )(+ — — I). Similarly in case (u) we 

find 1/4(1 ® (l+if7,)®(l + icrz))t TrA(|0)G(0|)(l ® 
(1 + icrz)® (1 -l-iCTz)) = (1 (g) Had®Had)V(l ® 
Had® Had). Hence the state obtained after tracing 
out 2 qubits is always separable, which implies that 
the relative entropy E^ vanishes. In addition, we have 
S{pa) — S{pb) and using the decomposition above 
we obtain S'(TrB(|0)G(0|)) = 2. The relative entropy 
of entanglement of the ensemble of pairs is given by 
Er{p'^) = Sfbcjes^bc as i?r = 1 for a fully entangled 
pair and Er = O for a separable state. We sum up the 
contributions of the (2) possible bipartitions to obtain 

9 mab<20M. (21) 

a,b-a<b 

Joining the two inequalities we get 

2 

M < -M. (22) 

and hence the procedure (fT5|) generates losses. On the 
other hand, breeding presented in Sec. IIVBI (or equiva- 
Icntly hashing) gives yield 1 for states of fidelity 1, and 
for a state of form ^ with / > 0.9877 we have F > 2/3. 
Hence, any 5-qubit ring-state arising from the application 
of global depolarizing noise and with fidelity / > 0.9877, 
can be purified more efficiently using the direct protocol. 

2. Minimal required fidelity 

In this section, we illustrate the process consisting in 
generating one entangled pair from a 5-qubit ring-state. 
We see that the minimal required fidelity is higher for 
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the bipartite strategy. We start with a state resulting 
from the apphcation of global white noise to the 5-qubit 
ring-state, given by Eq. [131 We define x = f — and 
rewrite the state as 

p = x\0)G{0\+^^^fil, (23) 

where G stands for the ring. To create a 2-qubit entan- 
gled pair from the initial state, we measure three consec- 
utive qubits in the eigenbasis of az ■ We remark that this 
is in fact an optimal strategy for states of the form Eq. 

when operating on a single copy. The resulting state 
of the remaining two qubits, after some local correction 
if the measurement result is 1, is given by 

p = x|0)G,(0| + ii^l, (24) 

where G2 is the graph consisting in two vertices con- 
nected by an edges. This state is equivalent up to lo- 
cal unitary operations to a Werner state with fidelity 
F — {3x+ l)/4. For x > 1/3 the state is distillable since 
bipartite entanglement purification protocols can be suc- 
cessfully applied. The state has positive partial transpose 
for a; < 1/3 which implies that it is not distillable. We 
thus have the condition x > 1/3 such that the state is 
distillable via a bipartite entanglement purification strat- 
egy. For the direct multiparty entanglement purification 
protocol proposed in this article, we (numerically) find a 
threshold x > 0.2 for states of the form Eq. Thus 
the minimum required fidelity for the multiparty strat- 
egy is significantly lower than for the bipartite strategy. 
Although here we illustrate the advantage of our method 
by the 5-qubit ring state, such advantages are expected 
for other graph states as well. 



where p is the reliability. As part of the purification pro- 
tocols, local one- and two-qubit unitary operations are 
employed which may be noisy. An imperfect operation is 
modeled by preceding the perfect operation {/^"^"^^ with 
the application of the noise superoperators £ from Eq. 
P5|l with parameter pi, i. e. the state is transformed as 

p ^ [/('^i'^^^) (4;i)4;'^)(p)) C/t('^i'^2), (26) 

We assume that the protocols are executed with the least 
possible number of operations to keep accumulated noise 
low. Hence, if a local two-qubit gate J/^j^"^^ is preceded 
by one-qubit gates [/}°^'' and 1/2'^^'^ we apply one com- 
bined unitary [/("i"^) = {/^^^'c/^^^'t/^a"'"-' which is sub- 
jected to noise only once. 

We compare a strategy using direct multiparty entan- 
glement purification (MEPP strategy) with a strategy us- 
ing bipartite entanglement purification (BEPP strategy) 
in the particular example of the 5-qubit ring-state, which 
is genuinely three-colorable and hence can be purified di- 
rectly only by means of the new protocol. To compare 
BEPP and MEPP strategies, we computed the maximal 
reachable fidelity Fmax, and the minimum required fi- 
delity Fmin- Ffnax IS the maximum value of fidelity to 
which a state of fidelity F > Fmin can be brought using 
the purification protocol. Note that F^ax is the same in 
a communication and in a static scenario. This comes 
from the fact that the maximal reachable fidelity does 
not depend on the initial state. That is, if the initial 
state is distillable, it is always possible to go to a state 
of fidehty Fmax- 



1. Communication scenario 



B. Performance under imperfect local operations 

After having shown the advantage of multipartite pu- 
rification with respect to bipartite purification in the 
static scenario, we turn to a more general setting includ- 
ing noisy local operations. We show here that when local 
operations are imperfect, multipartite purification can be 
advantageous also in the communication scenario. 

We model noise in the communication channels and in 
the local operations as follows. We study typical noise 
models, where the Kraus representation of the superop- 
erators is diagonal in the Pauli basis. This is a common 
and usually sufficiently general model [23 | (in particu- 
lar, any noisy channel can be brought to such a form by 
means of (probabilistic) local operations). We consider 
the depolarizing channel: 

4'')(p)=pp+i^(p + a('^)p4'^) + 

+ (25) 



Let us describe the MEPP and the BEPP strategies in 
the communication scenario. The initial states are dif- 
ferent in both strategies which renders the comparison 
non-obvious. We use the local noise equivalent (LNE), 
which is, for a state of fidelity F, the level of local de- 
polarizing noise defined as in Eq. (j25p . which has to be 
applied to the perfect state to obtain fidelity F. 

In the MEPP strategy, the setting is the following: 
5 parties A,B,C,D, and E are connected by depolariz- 
ing channels; the party A creates states locally and dis- 
tributes them to the four other parties. To use direct 
purification, in addition to the 5-qubit ring states, the 
party A needs to distribute 5-qubit cluster-states in three 
different ways, allowing the purification with respect to 
the three different colors. The purification is done in two 
steps. The parties first purify the three different clus- 
ter states up to their maximum reachable fidelity, before 
using them to purify the ring-states up to F^^^^ [pi), 
where pi gives the amount of local noise. The fact that 
F^^^^ [pi) is always reached for a state with F > 
^mm^^fe) is then used to compute F^^PP{pi). For 
a given amount of local noise pi, we vary the channel 
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FIG. 5: Local Noise Equivalent (LNE) corresponding to the 
maximal reachable fidelity and to the minimal required fi- 
delity as function of the amount of local noise pi for the 5- 
qubit ring state in the communication scenario. Given a state 
of fidelity F, the LNE is the level of local depolarizing noise 
which has to be applied to the perfect state to obtain fidelity 
F. The red solid line and the green dashed line stand for the 
LNE of the maximal reachable fidelity and the minimal re- 
quired fidelity respectively for the MEPP strategy. The same 
values are plotted for the BEPP strategy. The blue small 
dashed line is the LNE of a ring state teleported using 4 pu- 
rified pairs and the pink dotted line is th LNE corresponding 
to the minimal required fidelity of a pair obtained by send- 
ing one of its qubits through a depolarizing channel. Not that 
the maximal reachable fidelity is the same in a communication 
and a static scenario. We also plotted the maximal reachable 
fidelity for the alternative strategy where a 5 qubit cluster- 
state and a Bell pair are purified before being connected. This 
value is given by the light-blue dotted-dashed line. 



above which the 

MEPP ( 



noise q to find the threshold value qmin 
state can be purified up to F^J^f^ipi). F^{^^^{pi) is 
obtained by applying the depolarizing channel with noise 
parameter g„im to 4 of the qubits of the ring-state. 

In the BEPP strategy, the party A creates Bell pairs 
and distributes them to the other parties. The pairs 
are purified and then used to teleport locally created 5- 
qubit ring-states. We adopted a conservative scenario 
where one of the 5 parties creates the ring-states, de- 
creasing by this the number of Bell pairs needed to tele- 
port the states from 5 to 4. In addition, we assume that 
the teleportation process itself does not add additional 
imperfections/noise. Hence the actual value of F^^^^ 
is lower than our conservative estimate. To obtain the 
maximal reachable fidelity of multiparty entangled states 
P^aI^{Pi), we purify the Bell pairs up to F^^^^ipi) and 
use TV — 1 of them to teleport a locally created ring-state. 
The results presented in Fig. O show that in a commu- 
nication scenario, the minimum required fidelity as well 
as the threshold value pi , under which no purification is 
possible is always lower in the BEPP strategy. However, 
the maximal reachable fidelity is higher in the MEPP 
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FIG. 6: Maximal reachable fidelity and minimal required fi- 
delity as function of the amount of local noise for the 5-qubit 
ring state in the static scenario. The red solid line and the 
green dashed line stand for the maximal reachable fidelity 
and the minimal required fidelity respectively, for the MEPP 
strategy, while the blue small dashed line and the pink dotted 
line stand for the same quantities in the BEPP scenario. The 
value of Fmin plotted for the BEPP scenario is the minimum 
value of fidelity of the ring state, so that the pairs we extract 
from the ensemble of rings can be purified. 



scenario allowing us to get states of higher fidelity. 

We also studied the alternative strategy consisting in 
distributing two-colorable sub-graph states which are pu- 
rified and connected at the end. There are two different 
possibilities depending on the choice of the first color in 
the construction of the sub-graphs. If one chooses an 
ensemble of two qubits as first color, one gets a 5-qubit 
cluster-state and a 2-qubit cluster state (5-2 strategy), 
while if the first color contains only 1 qubit, one gets a 
3-qubits cluster-states plus a 4-qubit one (4-3 strategy). 
The classification of the different strategies by decreas- 
ing value of the maximal reachable fidelity is as follows: 
direct purification, 5-2 strategy, 4-3 strategy and finally 
the BEPP strategy. The order is inverted for the value 
of minimal required fidelity. 



2. Static scenario 

A similar calculation can be done for the static sce- 
nario. In this case, one is given an ensemble of distributed 

5-qubit ring-states described by p = nt=i ^p"^|0)g(0|- 
As noted in the previous subsection, the maximal reach- 
able fidelity is the same as in the communication scenario. 
In the MEPP strategy, the minimum required fidelity is 
similar in both scenarios. Indeed, the creation of the 3 
different cluster states can be done together with a first 
purification with respect to one of the colors. The be- 
haviour of the protocol in the static scenario therefore 
doesn't differ much from the one in the communication 



10 



scenario if one begins the purification of the two-colorable 
graph-state corresponding to graph gj, by sub-protocol 
Vj. However, the minimum required fidelity is differ- 
ent in the BEPP strategy. It is defined as the minimum 
value of fidelity of the ring-state, such that at least one 
distillable pair can be extracted from it. The results are 
presented in Fig.lHl The MEPP strategy clearly presents 
an advantage in terms of fidelity. The minimum required 
fidelity is lower and the maximal reachable fidelity larger 
for any value of reliability pi . 

VII. CONCLUSION 

In this paper, we have proposed multipartite entangle- 
ment purification (recurrence and breeding) protocols by 
which parties can distill arbitrary graph state directly. 
The work not only gives a complete systematic package 
for the construction of entanglement purification proto- 
cols for graph states, but also clarifies a special role of 
two-colorable graph states in reading out the parity in- 
formation of the stabilizer eigenvalues of any graph state. 
The latter property might open a new avenue to a "patch- 
work" purification of decohered qubits in resource entan- 



gled states for quantum information processing. We re- 
mark that very recently a similar idea, i.e. the usage of 
different shapes of graphs, has been utilized by Goyal, 
McCauley and Raussendorf to obtain purification proto- 
cols for two-colorable graph states with improved yield 
and scaling behavior [20|. 

We have considered two scenarios, namely (i) the static 
LOCC scenario and (ii) the communication scenario. Un- 
der ideal local operations, we have showed in the static 
scenario that our protocol gives a higher yield and a wider 
distillable regime (i.e., a smaller required fidelity for dis- 
tillability), compared with any bipartite strategy. Also, 
under noisy local operations, our protocol allows one to 
distill mixed states up to a higher achievable fidelity in 
both scenarios. 
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